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MAPPING THE GALACTIC CENTER WITH GRAVITATIONAL WAVE MEASUREMENTS USING PULSAR TIMING 
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ABSTRACT 

We examine the nHz gravitational wave (GW) foreground of stars and black holes (BHs) orbiting SgrA* in 
the Galactic Center. A cusp of stars and BHs generates a continuous GW spectrum below 40 nHz; individual 
BHs within 1 mpc to SgrA* stick out in the spectrum at higher GW frequencies. The GWs and gravitational 
near-field effects can be resolved by timing pulsars within a few pc of this region. Observations with the 
Square Kilometer Array (SKA) may be especially sensitive to intermediate mass black holes (IMBHs) in this 
region, if present. A 100ns-10/is timing accuracy is sufficient to detect BHs of mass 1000 M Q with pulsars at 
distance 0.1-1 pc in a 3 yr observation baseline. Unlike electromagnetic imaging techniques, the prospects for 
resolving individual objects through GW measurements improve closer to SgrA*, even if the number density 
of objects increases inwards steeply. Scattering by the interstellar medium will pose the biggest challenge for 
such observations. 

Subject headings: galaxies: nuclei - gravitational waves - pulsars 



1. INTRODUCTION 

There is a great ongoing effort to use pulsar timing arrays 
to detect gravitational waves (GWs) in the nHz frequency 
bands. GWs, if present, modify the exceptionally regular 
arrival times of pulses from radio pulsars. Observations of 
a correlated modulation in the time of arrivals (TOAs) of 
pulses from a network of highly stable milliseco nd pulsars 
(MSPs) across the sky can be us ed to detect GWs (iDetweilerl 
Ivan Haasteren et ail 1201 lb . Existing pulsar timing ar- 
rays (PTAs) utilize the brightest and most stable nearby MSPs 
in the Galaxy. 

At nHz frequencies, the GW background is expected to be 
dominate d by cosmological supe rmassi ye black hole (S MBH 
binaries (iRaiapopal & Romanilll995l: Uaffe & Backer! 120031 
IWvifhe & Loebl I2003t ISesana et al.1 120041) . Only a few 
studies considered the G W signal from nearby sources. 
lLommen & Backer! d2001l) showed that pulsar signals would 
be sensitive to a putative SMBH binary in the Galactic Cen- 
ter with mass ratio 0.06, however, such a binary would 
have other dynamical consequen ces which a r e not observed 
dYu&Tremaind 120031) . Further, Uenet et al.l (12004 showed 
that the nearby extragalactic source 3C 66B does not con- 
tain a massive SMBH binary, b ecause PTAs do not ob serve 
the expected GW mo dulation. iBlandford et all dl987l) and 
Ide Paolis et al.l (119961) examined if the variable gravitational 
field of nearby stars could be detected usi ng pulsar timing in 
the cores of globular clusters, and similarly. Uenet et all (120051) 
considered the possibility of detecting GWs from intermedi- 
ate mass black hole binary sources using pulsars in the clus- 
ter. 

In this paper, we examine the prospects for directly detect- 
ing the GW foreground and gravitational near-field effects 
of the Galactic Cent er (GC) using pulsars in the vicinity of 
this region (see also lRav & Kluzniakll 19941) . We estimate the 
foreground (in contrast with the cosmological background of 
GWs) generated by the dense population of stars and com- 
pact objects (COs) in the GC, including about 20,000 stellar 
mass black holes (BHs) dMorrislll993TlFreitag et al.ll2006al) 



and perhaps a few intermediate mass black holes (IMBHs) 
of mass 10 3 M Q dPortegies Zwart et aT] 120061) . As these ob- 
jects are much more massive than regular stars populating 
the GC, they segregate and settle to the core of the central 
star cluster. The emitted GW signal falls in the nHz range, 
well in the PTA frequency band. While this fore ground sig- 
nal ma y be faint at kpc distances in the Galaxy as iJenet et alj 
( 120041) discussed, it may exceed the GW background locally, 
near the GC. The present generation of PTAs, achieve an 
upper limit of the characteristic stochastic GW background 
amplitude of the order h c < 6 x 10" 15 at / ~ 1/ yr ~ 30nHz 
(Ivan Haasteren et aT1l201 11) . while the theoretical prediction is 
ft c ~9x 10" 16 (/yr)- 2 / 3 dKocsis & Sesan all20Tll and see §gj] 
below). We estimate the level of pulsar timing accuracy nec- 
essary (i) to constrain the mass of IMBHs in the GC using 
pulsars in the GC neighborhood to a level better than existing 
constraints, or (ii) to resolve the central cusp of stellar mass 
objects. 

While a large popula tion of pulsars is ex pected to reside 
in the Galactic Center dPfahl & Loebl 120041) . their detection 
is quite challenging. It requires high sensitivity (> 0.01 /iJy) 
at relatively high radio frequencies (> 10 GHz) where the 
pulse smearing due to t he scattering of the ISM is less severe 
(ILazio & Cordesl[T998h . S ince pulsars have a steep radio fre- 
quency spectrum, S u oc v~ X A to v 18 dKramer et al.|[l998l) . the 
high frequency of observation makes them very fa int and thus 
difficu lt to detect and time. Note however, that Kei th et aT] 
d201 ll) have successfully detected nine radio pulsars at a fre- 
quency of 17 GHz, including the detection of a millisecond 
pulsar and a magnetar with an indication that the spectral 
index may flatten above 10 GHz. The future Square Kilo- 
meter Array is expected to find se veral thousand regular pul 
sars, a nd a few MSPs in t he GC dCordes et alj|2004t ICordes 



2007[ ISmits et alj 120091: iMacquart et all 120 10t ISmits et al 



201 11) . High time resolution surveys rec ently found 5 mil - 



lisecond pulsars in mid galactic latitud es dBates et alj|201 11) 
3 within lOOpc of the Galactic Center d Johnston et alj 120061 : 
iDeneva et al.l2009l) . Based on the properties of nearby pulsars 
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and no ndetections in a targeted observation, iMacquart et alj 
( 120101) recently put an upper limit of 90 regular pulsars within 
the central 1 pc. It is likely that the timing accuracy of these 
pulsars in the GC will be much worse than those in the lo- 
cal Galactic neighborhood. However, the GW signal may be 
much stronger near the GCs to make a GW detection possible 
with a lower timing accuracy. 

We adopt geometrical units G = c = 1 , and suppress the G/c 2 
and G/c 3 factors to change mass to length or time units. 

2. BLACK HOLES IN THE GALACTIC CENTER 

2.1. Stellar mass BHs 

The number density of objects in a relaxed galactic cusp 
orbiting around an SMBH with semimajor axis a is 

3 — a f a \ a o 
n*(a)=— - AU— pc" 3 , (1) 

4-TT \pC J 

where AT* is the total number of objects within lpc, and 
a = 7 /4 for a stationary Bahcall- Wolf cusp dBahcall & Wolfl 
WWl iBinnev & Tremaind [2008). If the mass function in 
the cusp is dominated by heavy objects, the density of light 
and heavy objects relaxes to a profile with a ~ 3/2 and 
a ~ 7 /4, while in t he opposite case as steep as ~ 9/4 to 3 for 
the h eavy objects ( lAlexander & Hopma n 2009; Kesh et et"aT] 
120091) . These stationary profiles have a constant inward flux of 
objects which are eventually swallowed by the SMBH and re- 
plenished from the outside. We take the theoretically expected 
value N* = 20,000 f or stellar mass black holes of mass m± 
10 VI within lpc (iMorrisi 119931: iMiralda-Escude & Gould! 
12000b iFreitag et al.ll2006al) and assume a*,BH = 2. 

The eccentricity distribution for a relaxed thermal distribu- 
tion of an isotropic cusp is such that the number of objects in a 
de bi n is proportional to dN oc e de, independent of semimajor 
axis ( IBinnev & Tremainell2008l) . 

2.2. Intermediate mass BHs 

Intermediate mass black holes (IMBHs) are expected to 
be crea ted by the collapse o f Pop III stars in the early uni- 
verse dMadau & Reesll2~001l) . r unaway collisions of stars in 
the c o res of globular clust ers dPortegies Zwart & McMillan! 
l200l IFreitag etalJl2006bl) or th e mergers of stellar mass 
black holes dO'Learv et all 120061) . If globular clusters sink 
to the galactic nucleus due to dynamical friction, they are 
tidally stripped and deposit their IMBHs in the Galactic nu- 
cleus. Then the IMBHs settle to the i nner region of the nu- 
cleus by mass segregation with stars. iPortegies Zwart et alj 
(120061) predict that the inner lOpc of the GC hosts 50 IMBHs 
of mass M~ 10 3 M Q . 

There are very few unambiguous observations of in- 
termediate mass black holes (IMBHs) in the Universe 
(Miller & 'Coiber3l2004l) . Ultraluminous X -ray sources pro- 
vide the best observational candidates. In particular, HLX-1 
in ESO 243-49 is found to have a mass between 3 x 10 3 M Q < 
M < 3 x 10 5 M Q (iDavis et al.ll2011l) . In the Galactic Center, 
0.13pc projected distance from SgrA*, IRS 13E is a dense 
concentration of massive stars, which has been argued to host 
an IM BH of mass between 10 3 and 10 4 M© dMaillard et al.l 
120041) . however the observed acceleration constraints make 
an I MBH interpretat ion in IRS 13E presently unconvinc- 
ing dFritz et al.ll2010h . Astrometric observations of the ra- 
dio source SgrA* corresponding to the SMBH can be used to 
place an upper limit of the mass of an IMBH to M < 10 4 M Q 



in 5-500 mpc dReid & Brunfhalerll2004l) . Further, an IMBH 
coul d have served to deliver the observed young stars in the 
GC (lHansen & Milosavlievicl l200l iFuiii et "all 120091). eject 
hyper velocitv stars (lYu & Tremaina 12003b IGualandris et al. 
20051) . create a low-density core in the GC dBaumgardt et al. 
20061) . efficiently randomize th e eccentricity and o rientations 
of the observed S-star orbits (lMerritt et al.l |2009|). and may 
have contributed to the SMBH growth dPortegies Zwart et alj 
12006b . These dynamical arguments can be used to place in- 
dependent limits on t he existence and m ass of IMBHs in the 
Galactic Center (see iGenzel et al.l I2010L for a review). Fu- 
ture observation of the pericenter passage of the shortest pe- 
riod kno wn star S2 may improve this limit to a few x 10 3 M Q 
in 2018 dGualandris et al.l [20101) . and even better limits will 
be possible by imaging SgrA* with Event Horizon Telescope 
(EHT), a milUm eter/submillimeter very long baseline interfer- 
ometer (VLBI) dBroderick et alJl2011h . We examine whether 
pulsar timing could detect an IMBH with parameters not ex- 
cluded by existing observations, or be used to place indepen- 
dent limits. 

2.3. Loss cone 

A depleted region is formed in phase space if objects are 
removed at a rate faster than they are replenished from outside 
by inward diffusion. In this region, Eq. (Q3 is no longer valid. 
The dominant source of removing stars is tidal disruption or 
physical collisions, while for BHs it is GW capture by the 
SMBH. 

The objects with initial conditions (a,e) fall in and merge 
with the SMBH due to GW emission in a time 

s 4 

U g = ^ ; ^(l-e 2 ) 7/2 = 36/ t mJ 1 a 4 npc (l- e 2 ) 7 / 2 Gyr. (2) 

where a mpc = a/mpc, m^ = to*/10 3 Mq, an d 1 < k < 1. 8 is a 
weakly dependent function of eccentricity dPetersll 19641) . As- 
suming e ^> and < f mg , the minimum semimajor axis is 

a lc = 0.4 mpc x in^d-e)- 1 ^ , (3) 

where the inward diffusion time is parameterized as fd,9 = 
fd/10 9 yr, and we assumed / ~ f p ~ (1 - e)~ 3//2 /orb ( see Eq. 
( IA5b ). For stellar mass BHs, the inward diffusion time is re- 
lated to two-body relaxation dBinney & Trem aine 200 8]). De- 
pendin g on the number and masses of BHs. lO'Learv et alj 
d2009l) find that 0.1 < f d , 9 < 10. For IMBHs, the inward dif- 
fusion is due to dynamical friction and the scattering of stars. 
This process is initially faster than the relaxation timescale, 
but then slows down (t^g ~ 10) as stars on crossing orbits 
are ejected by the IMBH dGualandris & Merrittl 120091) . The 
number density inside a\ c is expected to be much less than 
that of Eq. ([T). In such a state, the eccentricity of an 
IMBH is increased dMatsubavashi et alj l2007t ISesanal 120101: 
llwasawa et alJl20TlT) . However, a possible triaxiality of the 
cluster might result in the refilli ng of the loss-cone and shorter 
inward migration timescale s (Khan et al] 1201 lb iPreto et aT] 
l20llt IGualandris & Merrittll2012l) . 

3. GRAVITATIONAL WAVES FROM THE GALACTIC CENTER 

We start by reviewing the essential formulas to derive the 
GWs generated by a population of binaries with circular or- 
bits, then turn to the general eccentric case. We discuss other 
details of the spectrum, regarding the high frequency cutoff 
and splitting into discrete peaks, at the end of the section. 
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3.1. Unresolved circular sources 

The GW frequency for a circular orbit is twice the orbital 
frequency, / = 2/ OT b, and the corresponding orbital radius is 



r(f)=M.(nM.f) 



-2/3 . 



2.6 / 8 2/3 mpc. 



(4) 



where = //(10 8 Hz). In the last equality, the mass of 
the central SMBH in SgrA* is taken as M. = 4.3 x 10 6 M Q 
dGillessen et al .1120091) . To put in context, 20 mpc (i.e. about 
4,100 AU) is the distance to the observed S-stars, the in- 
nermost star, S2, has semimajor axis 4mpc and pericenter 
0.8 mpc, and the stellar disk of mas sive young stars extends 
down to 30mpc dGenzel et al.ll2010l) . 

The RMS strain generated by an object of mass orbiting 
around a SMBH of mass M, on a circular orbit at distance D 
from the source in one GW cycle is 



32 M.m* 

5~ Dr(f) 



: 8.8 x 10- 15 m 3 £>pc/f ; 



2/3 



(5) 



where D pc = D /pc, mi = m* /(10 3 M Q ), and the index wiU 

stand for zero eccentricity. The y/32/5 prefactor accounts for 
RMS averaging the GW strain over orientation. 

The GW strain of many independent sources with the same 
frequency but random phase adds quadratically. For a sig- 
nal observed for time T, the spectral resolution is Af = l/T. 
Therefore the number of sources with overlapping frequencies 
is 



A »r dN 

AN= — 
dr 



dr 
df 



The number of GW cycles observed in time T is fT. 
integrated GW signal with frequency / is 



h 2 c = AN(fT)hl = 



dN 
dr 



dr_ 
df 



m= 



dN 



hi 



(6) 



The 



(7) 



also called characteristic spectral amplitude (lPhinney||2001l) . 
Therefore the characteristic GW amplitude from a population 
of sources can be interpreted as the RMS GW strain in a log- 
arithmic frequency bin. Here dN /dr is the number of sources 
in a spherical shell and dr/df is to be evaluated using Eq. @. 
Note that the final formula is independent of rQ 

We can express h c with the number density of objects using 
dN /dr = Airr 2 n+(r), as 

Afr) = hi = ^ ^ (.M.f)^ n Mf)] ■ (8) 

Eq. ^ gives the root-mean-square (RMS) GW signal level 
drawing the stars or BHs from a density profile n*(r). 
The equation shows that the spectral shape is different 
from f~ 2 / 3 , de scribing the RMS cosmological background 
(lPhinneyl2001h . Note that h c is proportional to the RMS mass 
of objects which may exceed the mean for a multimass popu- 
lation. The GW signal for any one realization is well approx- 
imated by the RMS if AN > 1 . Close to the center (which 
corresponds to high /), AN is small, and the GW spectrum 
becomes spiky, which we discuss in § l3.3l below. 

The GW foreground may be different from the above esti- 
mate for sources with significant eccentricity, for sources on 

1 This is true as long as the source distribution is effectively continuous 
in frequency space, so that indiv idua l sour ces a re unresolved, but not for 
resolved discrete sources (see Eq. OD and jj-FOlbelow). 



unbound orbits, or if they are non-periodic or evolve secularly 
during the observation, or if the population is anisotropic. We 
elaborate on the corresponding effects on the GWs in the sub- 
sections below, and discuss the characteristic GW frequency 
where the spectrum would be affected by the loss cone and 
small number statistics. 

3.2. Eccentric periodic and burst sources 

Here we summarize the modification of the above estimates 
due to eccentricity, and refer the reader to Appendix for 
details. 

For eccentric sources, the GW spectrum of individual 
sources is no longer peaked around / = 2/ or b. For mildly 
eccentric sources (e < 0.3), the GW spectrum is spiky, with 
discrete upper harmonics, / = nf or \,, with decreasing ampli- 
tude for n > 2. The observed width of each spectral peak is 
Af = l/T. For larger eccentricities, the upper harmonics are 
stronger than the n = 2 mode, and the peak frequency corre- 
sponds to the inverse pericenter passage timescale f p , where 
90% of the GW power is between 0.2 f p < f < 3f p . 

The GWs are in the PTA frequency band if the pericen- 
ter passage timescale is less than the observation time, T. 
We distinguish two types of sources: periodic and GW burst 
sources^ Periodic sources have orbital periods shorter than 
T, or semimajor axis 



^per — 5 .7 7^jq 



/3 



mpc . 



(9) 



At relatively low frequencies where small number and loss 
cone effects do not kick in, these sources generate a continu- 
ous GW spectrum which is similar to that given by the circu- 
lar formula Eq. $13[ within a factor 2, ass umin g an i sotro pic 
thermal eccentricity distribution (see Eq. ( IA12b and (I A 1 3t in 
AppendixlAl. 

Sources with semimajor axes larger than Eq. (O generate a 
GW burst during pericenter passage in the PTA band if their 
pericenter distance at close approach is less than Eq. (0. Only 
a small fraction of these burst sources contribute to the PTA 
measurements, the ones which are within time T from peri- 
center passage along their orbits during the observation. In 
AppendixlAl we show that the stochastic background of burst 
sources is expected to be typically much weaker than the peri- 
odic sources, unless there are IMBHs on wide eccentric orbits. 

3.3. Individually resolvable sources 

The GW foreground generated by a population of objects is 
smooth if the average number per Af frequency bin satisfies 
(AN) 3> 1. If the number density follows Eq. (Q~|), and the 
total number within lpc is normalized as iV* = A^/(2 x 10 4 ), 
then the spectrum becomes spiky ((AN) < 1) above 



/ res =4.2xl0- 8 Hzxl0 9(Q - 2) / (9 - 2Q) 



(3-a)2V* 



3/(9-2a) 



(10) 

where 7\o = T/lOyr, and we used Eqs. (HJi and © for circular 
orbits. The corresponding orbital radius is 



1.0mpcxl0 6(a - 2)/(2a - 9) 



1 2/(9-2a) 



(3-a)JV* 



(11) 



2 There are also "repeated burst" sources whi ch are highly eccentric with 
orbital periods less than T, which satisfy Eq. (9) I Kocsis & Levin 201 1). We 
do not distinguish repeated bursts from periodic sources here. 
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Sources within r res generate distinct spectral peaks above fre- 
quency / res . We refer to these sources as resolvable. 
The total number of resolvable sources is 



N m = 



n*0-)47rr 2 dr = jV* ( — 
pc 



(12) 



: 20 x io 9 («- 2 >A 9 - 2 «>jv* 



(3-aW* 



2(3-a)/(9-2ct) 



Note that / res , r res , and N res are exponentially sensitive to the 
density exponent: in a 10 year observation of N+ = 20,000 
BHs, / res = (19,42, 110)nHz, r res = (1.7, 1.0,0.53)mpc, 
and = (7,20,70) for a = (1.75,2,2.25), respectively. 
Similarly, for N+ = 50 IMBHs, / res = (2.0,3.8, 6.4)nHz, 
'"res = (7.6,5.1,3.5)mpc, and N res = (1,3,12) for a = 
(2.25,2.5,2.75), respectively. If the supply of objects by two 
body relaxati on i s very slow, t j 3> 1 Gyr, and the loss cone is 
empty (see § 12.31 ), then the number of resolvable sources can 
be much less. The GW spectrum transitions from continuous 
to discrete inside the PTA frequency band, and the expected 
number of resolvable sources is typically non-negligible. 

For / > / les , the GW spectrum is h c {f) = 0, except for 
distinct frequency bins which include one resolvable source 
each. For the latter, the net GW signal amplitude in time T is 

KAf) = (fT^hoif) = 1-6 x l<r H m 3 D£Ttf 2 fl /6 . (13) 

Therefore h c \{f) increases quickly for individual sources 
with decreasing orbital radius or increasing frequency. The 
prospects for detecting individual sources closer than r res to 
SgrA* improves because the GW amplitude increases, while 
the confusion noise per frequency bin decreases (see Eq. ( l23l 
below for the co rresp onding timing residual). 

Note that Eq. (TT~3b corresponds to the GW polarization av- 
eraged integrated signaQ for stationary circular orbits. We 
derive h c {f) for eccentric periodic and burst sources in Ap- 
pendix [A] In that case, the GW power of each source is dis- 
tributed over many harmonics, making the detection of indi- 
vidual sources potentially more challenging. However, GR 
pericenter precession leads to an overall amplitude modu- 
lation of the relative weights of th e two GW polarizations 
dKonigsdorffer & Gopakumarl2006l) . This effect increases in- 
ward, which may improve the prospects for separating indi- 
vidual sources. 

3.4. Maximum GW frequency 

What is the maximum GW frequency for objects outside of 
the loss-cone, for which the GW merger is still longer than the 
relaxation timescale? The maximal emitted GW frequency 
corresponding to Eq. (01 is 



/ic = 1.6 x 10"' Hz x m 3 3/8 (l-e)- J/16 f, 



-3/16,-3/8 



d.9 



(14) 



where note that the dimensionless diffusion time is 0.1 < 
fd.9 < 10 for multimass two-body relaxation models (see 
§ I2.31 l. and we assumed / ~ / p ~ (1 -e)" 3 ' 2 / or b (see Eq. 
( IA5I )). The GW spectrum is expected to terminate above f c . 
Eq. (TBI) implies that 1OM BHs generate a GW background 
with a maximum frequency of around 10~ 6 Hz = (12 day)" 1 
for mildly eccentric sources, and somewhat larger for more 

3 More precisely, for circular sources, — (J^ (h\+h 2 x )fdl), where (.) 
corresponds to averaging over source inclination. Note that throughout this 
manuscript h c (f) is in dimensionless units, counts per frequency bin. 



eccentric orbits. For IMBHs with m = 1O 3 M , t d ~ lOGyr, 
and e = 0.9, the corresponding cutoff frequency is around 
2 x 10~ 7 Hz = (60day) _1 . Sources outside the loss cone are 
in the PTA GW frequency band. 

3.5. Frequency evolution 

The above treatment is valid only if the frequency shift is 
negligible during the observation, which we discuss next. 

The timescale on which the GW frequency evolves, is 
/// ~ td, where f<j is the timescale on which sources drift in- 
wards. The frequency drift is then fT ~ fT jt&. Compared to 
the frequency resolution Af = l/T during the observation, 



fT_ fT_ 



3x lO-VsTfo^- 



(15) 



The frequency drift is negligible due to GW emission or two- 
body relaxation for sources outside the loss cone. 

4. GRAVITATIONAL WAVES AND PULSAR TIMING 

The cosmological GW background and its detectability 
with PTAs can be summarized as follows. 

4.1. Stochastic cosmological background 

The cosmological GW background constitutes an astro- 
physical noise for measuring the GWs of objects orbit- 
ing SgrA* in the GC. At these frequencies the stochastic 
GW background is expected to be dominated by SMBH 
binary inspirals, with a characteristic amplitude roughly 
(lKocsis&Sesanall2011l) 



hc. = l. 



xl0- 15 / 8 " 2/3 . 



(16) 



The GW background may be suppressed below the stochas- 
tic level of Eq. (THjil by a factor 2-4 due to gas effects at fy < 1 
dKocsis & S esana 20 1H), and bec ause of small number statis- 
tics at fs > 1 dSesana et alj|2008l) . On the other hand, individ- 
ual cosmological SMBH binary sources may stick out above 
this background by cha nce if they happen to be very massive 
and relatively close-by dSesan a et al. 20091). 

4.2. Sensitivity level of PTAs 

If pulsars are observed repeatedly in At time intervals for 
a total time T, the spectral resolution is A/ = l/T, and the 
observable frequency range is l/T < f< 2/ At . For T = lOyr, 
At = 1 week, therefore 3 x 10" 9 Hz < / < 3 x 10" 6 Hz. 

To derive the angular sensitivity of PTAs, let the unit binary 
orbit normal, the vectors pointing from the Earth toward the 
pulsar, and from the binary to the pulsar be respectively, L a = 
(sin l cos A, sin i sin A , cos c), p a = (sin p cos ip, sin p sin ip , cos /j.) 
and ¥' = (0,0, 1). In the far-field radiation zone of a circular 
binary, the metric perturbation is h a h = h* b {f)+h* b (f), where 



h + ab (t) = 2[l+(L%) 2 ]Asm<f>(t)e > a 
h* b (t) = 4L c k c Acos<5>{t)e + ab , 



(17) 
(18) 



where kb gives the GW propagation direction, $(?) = 2ivft is 
the GW phase, and A = M m m+/(D r) is the GW strain ampli- 
tude without constant prefactors, which we assume is much 
larger at the pulsar than at the Earth. Note that D is the 
distance between M. and the pulsar, while r is the distance 
between M. and m*. We orient the coordinate system so 
*at e* b = e x a e y b + e x b el and e + ab = with «J = (1,0,0), 



5 



e y a = (0, 1,0). Assuming that the pulsar is much closer to the 
source than the Earth, the TOAs of pulses change as 



5t = z(t)dt, 



(19) 



where z(t) = [f(t)— fol/^o = z+(t)+z x (t), where v(f) is the ob- 
served pulsar frequency, Uq is its unperturbed value (i.e. ne- 
glecting GWs), and 



1 p a p b 

h* b = (1 - cos /x) cos 2ip (1 + cos 2 t)-<4 cos $ 



Zx = 



2 1 + /) fl /t fl 
1 p a p b 



(20) 



/z^, = (1 -cos^)sin2?/>(2cosi)„4sin$ (21) 

Z, i ~\~ J} Kfl 

dDetweilerl 1979b lFinnl2009t lCornishl200"9T) . Note that unlike 
for cosmological sources the Earth-term is not included be- 
cause it is negligible in the timing residual for nearby pulsars. 
Averaging the redshift factor over an isotropic distribution of 
p a and L a , and over one GW cycle, 



ZGW= \/(Z+) + (Zx) = 



M * m * 1 ; tf\ m\ 
— = —^nn(f), (22) 

15 Dr(f) V6 



where h (f) is the RMS GW strain defined in Eq. ©. The 
variations in the TOAs are given by the integrated change, 
Zgw /2nf. For an observation time T, the sensitivity increases 
with the number of observed cycles as (/T) 1 /2 , 



= 100m 3 D- ( !r i 1 /2 / 8 1/6 ns. 



2irf J v^27r/ 

(23) 

Thus, instead of the integrated strain h / (2irf), the timing SNR 
is proportional to h c /(2-irf), which incorporates the (/T) 1/2 
factor. In the last equality, we assumed the contribution of a 
single source only, Eq. ( fl3l ), but in general h c can represent 
many sources per bin, given by Eq. ((SJ. Conversely, for fixed 
St timing accuracy, the detectable characteristic GW ampli- 
tude for an average binary orientation and pulsar direction, is 

/j,,pta(/) = V6(2nf)St = 1.1 x lO" 17 ^ , (24) 

independent of T, where St\o = <5f/10ns. The sensitivity is 
better in certain orientations: by an average factor s/3 if the 
binary-pulsar-Earth angle is 90° (with random binary inclina- 
tion), by y/5/2 if the pulsar is on the binary orbital axis (with 
random pulsar-Earth direction), and by y/3 x y/5/2 = 2.7 if 
both conditions hold. 

As of 2009, 13 pulsars were known with timing accu- 
racy better than St ~ 250 ns, the best ones reaching 20-30 ns 
dHobbs et al.l 120 id) . The huge collecting area of the future 
Square Kilometer Array (SKA) implies that it will be able to 
time many hundreds of pulsars at the level currently achieved 
for a few pulsars, and w ill perhaps find on es with much better 
characteristics (see also lJenet et alj|2009l) . However, the tim- 
ing accuracy may be much worse for fainter sources further 
away (see §|6]below). 

4.3. Near-field effects 

Up to this point, we assumed that the pulsar signal is modu- 
lated by the leading order radiative GW terms, and neglected 
other post-Newtonian near-field effects. This is justified if the 



binary to pulsar distance is D p 3> A, where 



1 



A= — = 0.16/ 8 - 1 pc. 
Zirf 



(25) 



is the reduced GW wavelength. Thus, near-field effects will 
be negligible for /g > 1 for D p > 1 pc, but they are significant 
at lower frequencies for pulsars closer to the center. 

Let us estimate the order-of-magnitude of the leading or- 
der near-field corrections for the particular geometry of this 
problem. The variations in the TOAs are approximately de- 
termined by the metric g M „ at the pulsar. The near-zone met- 
ric induces a variation in the pulsar po sition and velocit y, and 
modulates the propagation of pulses dJenet et al.l 120051) . We 
resort to simple order-of-magni tude estimates neglecting pos - 
sibly important lensing effects dWex et alJll996uFinnll2009t) . 
In the post-Newtonian (PN) approximation, the metric is prac- 
tically a power series in m,-/r,-, ri/rj, and v,, where m, are 
the masses (i.e. M, or mj, and rj ar e distance parameters 
(i.e. r or D), and v, are the velocities (|Blanchetet al.l 119981 : 
lAlvill2000t iJohnson-McDaniel et al.ll2009h . The constant co- 
efficients in this power series are typically order 1 . We look 
for periodically varying terms in the 2PN metric of the bi- 
nary which might be comparable to /jo(/) at the pulsar, and 
assume the corresponding correction to the gravitational red- 
shift and Doppler shift, collectively called Einstein delay, is 
proportional to these terms. 

For this estimate we restrict to circular binary orbits, where 
m* <C M, <C r <C D, with m* ~ 10 3 M o , r ~ mpc, D ~ pc. In 
this regime, 



m* M. M. , 2/3 

— <C <C — -~7.3x 10 / 8 ' . 
r D r 



•5xl0- 3 / 8 1/3 : 



(26) 



We are not interested in terms that are not time-varying such 
as M,/D or M;/D 2 , since these generate a constant gravita- 
tional redshift. Remarkably, the mass dipole terms m^r^/D 2 
andM.r./D 2 are much larger than the standard radiative term 
m±M,/rD at these distances, but these terms cancel out in the 
center of mass frame and are not present in the metric. (Here 
r. =M.r/(M. + ra i ) and = m*r/(M.+m*) give the vectors 
to the binary components relative to the barycenter, r is the 
binary separation, D is a field point.) However, the pulsar per- 
ceives the binary components at their retarded positions and 
so the dipole terms do not cancel out exactly unless the bi- 
nary line-of-sight velocity is zero. The current dipole terms 
m*(v* • ^Xr* ■ k)D~ 2 are indeed present in the goo, go,, and gu 
components of t he binary metric, (see Eq. (2. 15 ) in lAlvil2000l : 
and Eq. (6.4) in lJohnson-McDaniel et"alll2009l) . This leads to 
an orientation-averaged TOA correction 



m+vr y/JT 2 -5/6™i/2 

St cd —^^r = 30m 3 D ^f s J, n ' ns. 



D 2 2tt/ —"•^vo-'t 'io 
The next correction is the mass quadrupole, 



St 



tnq 



D 3 2irf 



n r->-3 r-1 l/6~l/2 

9m 3 £) pc / 8 ' T 10 ' ns 



(27) 



(28) 



Now consider terms higher order in the masses. The 
leading order term here is the quadrupolar radiation term 
M,m±/(Dr) ~ m+v 2 /D. which leads to TOA residuals i5?gw 
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FIG. 1 . — The characteristic strain amplitude h c vs. GW frequency / for a pulsar at 1 pc (left) and 0.1 pc (right) from the Galactic Center. Dotted black lines 
show the orientation-averaged h c for fixed timing residuals. Blue and red dashed lines show respectively the orientation-averaged timing residuals caused by 
individual stellar BHs and IMBHs on circular orbits in a 10 yr observation, assuming 1 object per frequency bin. Magenta lines show the timing residuals for a 
random realization of 1OM0 stellar BHs in the cluster (20,000 BHs within 1 pc with number density tx r~ 2 ). Typically only few bins are occupied at high /, 
generating a spiky signal with distinct sources. At lower /, many sources overlap to create a continuous spectrum. The RMS average over many realization of 
the cluster is shown by solid blue. The cosmological RMS stochastic GW background (green dashed) can be much smaller. 

of Eq. d23l i. All other terms in the metric are smaller by posi- 
tive integer powers of the small parameters in Eq. (l26l l. 

In addition to the Einstein delay discussed above, the tidal 
effects of the binary also induces epicyclic motion in the pul- 
sar, which leads to variations in the path length to Earth. 
This leads to a timing variation, St^, analogous to the Roe- 
mer delay in pulsar binaries relative to the binary barycenter. 
We estimate r p , the radius of the pulsar's epicyclic motion 
around the guiding center of its mean motion. We assume 
that the pulsar exhibits oscillations on the forcing period, so 
that its angular velocity matches the angular frequency of 
the forcing, tu. The centripetal acceleration is then r p uj 2 . 
We consider the case where the forcing is due to a current 
dipole or a mass quadrupole, |/^d,R | / m p ~ »?*vrsini/D 3 and 
|-Fmq,R|/'%> ~ »?*r 2 sini/D 4 . The time lag is St& ~ 2r p cos/x in 
a half cycle, and we assume the measurement improves with 
ifT) 1 ' 2 for longer observations. Neglecting factors of order 
unity, we get 



and the net timing residual can be estimated by scaling with 
^d/V/dln/ using Eq. ©. 

5. RESULTS 

We can now combine the results above to draw conclusions 
on the detectability of GWs from the GC with pulsars in its 
neighborhood. From Eq. $23[ , the distance within which a 
PTA could measure the GWs of an individual source with a 
fixed timing precision St = 10<5?ions is 

Ds^Um.St^ff^pc. (31) 

Eqs. ( fT3] ) and ( fl6l l show that the GWs from an individual BH 
in the GCs rises above the stochastic GW background within 
a distance 



o _ T l/2 Jl/6 

8.7OT 3 r i0 ' / g pc. 



(32) 



St, 



cd,R 



St 



D 3 



m+r" 



mq.R 




: 5 m^D 



pc 



3 ,-11/6 1/2 
J$ 1 10 nb ' 



D 4 (2Trff 



(29) 



(30) 



These estimates are consistent with that of Uenet et al.l ( 12005b . 
They found timing residuals corresponding to <5f mq .R can be 
500 ns for a pulsar inside a globular cluster with D pc = 0.03, 
/g = 0.3, 7io = 1 for a 10M o BH orbiting around an 10 3 M Q 
IMBH. The residuals are much larger in the near-field zone of 
the Galactic nucleus where the central object is a SMBH. 

In conclusion, the Einstein delay terms <5? cc j and St mq given 
by Eqs. (|27h28T i are much larger than the Roemer delay <5f cc j,R 



and <5f mq R for D > pc and / > lOnHz. In particular, <5? c( j is 
comparable to the standard radiative GW modulation, <5?gw at 
lOnHz. The gravitational near field timing residuals of indi- 
vidual sources become much larger for pulsars at D ~ O.lpc, 
at the level St c d ~ St mq ^ ~ lOm^fis at / ~ 10m3nHz, and 
even larger for smaller /. For stellar mass BHs, the individ- 
ual contributions of residuals are overlapping in Fourier space 
where near-field effects are most important (see Eq. ([Toll). 



A pulsar within D§, and £>b g to the GC could be used to 
detect GWs from individual objects in the GC. These esti- 
mates are conservative. First, they assume an average or- 
bital and pulsar orientation, D\, and Ds, might be 2-3 x larger 
for certain orientations. Second, Eq. d32b assumes that small 
number statistics and gas effects are negligible in cosmologi- 
cal SMBH mergers, and might overestima te the background 
dSesana et alJl2008l:lKocsis & Sesanall20ll FI 

Figure [T] shows the orientation averaged characteristic GW 
amplitude for a 10 year observation, incorporating the addi- 
tional near-field effects using Eqs. d24l i and ( 1271430b . The 
lower and top .Y-axis shows the GW frequency and orbital 
radius for circular sources. The magenta curve displays the 
spectrum of timing residuals for a Monte Carlo realization of 
a population of 20,000 BHs within 1 pc with number density 
n(r) ex r~ 2 and random orientation with mass 1OM on cir- 
cular orbits. The solid blue line shows the RMS foreground 
of a cusp of stellar mass BH averaged over different realiza- 
tions. The spectrum separates into distinct spectral spikes at 
higher frequencies with RMS maxima shown by the dashed 
line. Figure [2] shows the GW amplitude of a realization 

4 Note that the signal from individual cosmological SMBHs can also be 
enhanced by a factor 2-3 for certain orientations, but this effect is washed out 
for the unresolved stochastic background. 
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FIG. 2. — Same as Fig.[T]but for 1000 BHs of mass 30Mq within 1 pc, with a steeper density profile r 2 ' 5 . There are more resolvable sources in this case. 



of 1000 BHs with mass 30 Mg in lpc with a steeper num- 
ber density profile n(r) oc r~ 2 5 . Despite of the smaller overall 
number of sources in the cluster, there are many more resolv- 
able sources within 1 mpc in this case. 

The net background level in Figures[TJand[2]is conservative, 
as the background is sensitive to the RMS mass of objects in 
the cluster, which may significantly exceed the mean. Fur- 
thermore, individual sources may generate up to ^ 3 x larger 
timing residuals for certain orientations. Timing pulsars at 
a distance 0.1-lpc from the Galactic Center with 100ns- 
10 fis timing precision can be used to detect IMBHs of mass 
10 3 M Q , if they exist within 6 mpc of SgrA*. A popula- 
tion of stellar BHs within 2-5 mpc generates timing variations 
greater than 100ns-10/^s for a pulsar within O.lpc. However 
it will be much more difficult to individually resolve 10M Q 
stellar BHs within 1 mpc to SgrA*, which requires an extreme 
1-20 ns timing precision. 

Figure [3] shows the characteristic GW spectra for a Monte 
Carlo realization of an eccentric population of 10M Q BHs 
with an isotropic thermal eccentricity distribution. Magenta 
and cyan lines show the contribution of periodic and burst 
sources, respectively. The dashed lines show h c (f) for indi- 
vidual circular sources, and the solid blue lines represent the 
RMS GW level for a population of circular sources for com- 
parison. The net GW spectrum of burst sources is typically 
less than the level of periodic sources. The figure verifies the 
analytical calculations of Appendix [A] the continuous low- 
frequency spectrum of periodic sources is indeed comparable 
to the circular level, modulo a weakly frequency dependent 
constant between 0.4 < K < 0.8. Note that for clarity, we are 
not including the gravitational near field effects here, which 
would dominate over the continuous low-frequency GW spec- 
trum for D ~ 0. 1 pc as shown in Figure [TJ 

The transition to a spiky spectra happens at somewhat larger 
frequencies for eccentric sources (c.f. Figs. [TJ and |3). Since 
individual sources generate GW spectra with many orbital 
upper harmonics, the net high frequency spectrum is more 
complicated than in the circular case. The orientation av- 
eraged spectral level for individual sources is typically less 
than St < 200nsf7Z3Dp C 1 per frequency bin, however, the to- 
tal root-sum-squared signal of all upper harmonics of indi- 
vidual sources is much larger than the level of a single bin. 
In Appendix lAl we show that the SNR for detecting individ- 
ual sources with a matched filter is comparable for eccentric 



and circular sources. In this sense, the dashed lines in Fig. [3] 
are representative of the total timing residual of individual 
sources as a function of pericenter frequency for arbitrary ec- 
centricity. However, the full spectrum is rich in narrow fea- 
tures for eccentric sources and pericenter precession slowly 
modulates the amplitude of the timing residuals in individual 
pulsars. Both of these features could help to separate indi- 
vidual eccentric sources from the signal of other sources in 
GC. 

6. DISCUSSION 

We have shown that pulsars within a few pc distance from 
the GC offer a unique probe to identify IMBHs and stellar 
BHs orbiting around SgrA*. An IMBH, if present in the GC, 
sinks to orbital radii corresponding to the pulsar timing fre- 
quency bands. Depending on the binary orientation, orbital 
frequency, and pulsar distance, the GWs and gravitational 
near- field effects modulate the TOAs by a few ns to 100/is 
for these sources with masses between 10M Q and 10 4 M Q . 
Based on the GW spectral features, the signals of more than 
one IMBH (up to tens, if present) could be individually re- 
solved and isolated from the fainter signal of stars and stellar 
mass objects and the cosmological GW background. 

This observational probe is complimentary to EM mea- 
surements with different systematics. GWs are generated 
by all gravitating objects, including those that are black and 
undetectable in EM bands. GWs escape the galactic nu- 
cleus without any dissipation or dispersion. Closer to SgrA*, 
sources generate stronger, higher frequency GWs, the num- 
ber of observable cycles increases, and the number of objects 
per frequency bin decreases. Thus, unlike EM imaging tech- 
niques, the prospects for detecting and resolving individual 
objects through GW measurements improve closer in towards 
SgrA*, even if the number density of objects increases in- 
wards steeply. Furthermore, the gravitational effects are pro- 
portional to the RMS mass of objects, making the measure- 
ment more sensitive to individual higher mass objects in the 
distribution, even if the total mass of lighter objects is some- 
what larger on comparable radius orbits. 

Repeated pulsar timing observations over a few year base- 
line could reveal a detailed census of BHs in the inner mpc 
of the GC. We have shown that the GW foreground of stel- 
lar BHs rises above the cosmological background and sep- 
arates into distinct peaks above 40nHz, corresponding to a 
GW period of less than 1 yr or orbital separations less than 
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FIG. 3. — Characteristic spectral amplitude h c (f ) for timing a pulsar at 1 pc (left) and 0.1 pc (right) from the GC for a population of IOMq BHs with isotropic 
thermal eccentricity distribution. The contribution of periodic sources are shown with magenta, burst sources with cyan. For comparison, the dashed blue and red 
lines show the RMS GW level for individual circular sources, and solid blue lines show the RMS level from a population of circular 1OM0 sources (see sections 
2.1 and 2.2). Note that timing variations due to gravitational near field effects are not shown, which dominate within 3 X 10~ 8 Hz (see Figures [T]and|2j. 



1 mpc. Based on a simple estimate using circular orbits, we 
found that the total number of individually resolvable BHs is 
between 7-70, depending on the radial number density dis- 
tribution exponent r~ a between r" 1 75 and r~ 2 25 . These ob- 
servations are therefore exponentially sensitive to a, capable 
of tes ting the theory of strong mass segregation ( Keshet et al. 
12001 . 

Eccentricity complicates the spectral shape of resolvable 
sources by adding upper frequency harmonics. Although the 
timing residuals in individual frequency bins is suppressed by 
this effect relative to circular orbits, the total timing residual 
SNR with a matched filter is comparable for eccentric sources. 
At lower frequencies, a population of objects on eccentric or- 
bits generates a stochastic GW foreground with a similar spec- 
tral shape and a comparable amplitude, as a circular popula- 
tion. Objects on larger-semimajor axis, eccentric orbits gen- 
erate GW bursts during pericenter passage near SgrA*. The 
stochastic GW burst signal is much less than the level of peri- 
odic sources. 

This analysis hinges on the assumption that future surveys 
will discover pulsars near the GC that can be timed to the suf- 
ficient accuracy. Most of the observed S-stars within a few 
mpc and the young O/B stars in the G C will eventually tur n 
into pulsars in a supernova explosion (I'I'ahl & Loebl 120041) . 
Based on the age and number of these stars, and assuming 
that we don't live in a special time, one might expect more 
than 10 4 pulsars in the GC. Some of these may become MSPs, 
and may be beamed to wards us to be dete ctable with future 
SKA-type instruments dCordes et al.l 12004b . They might be 
expected to segregate to the outskirts o f the GC on a Gyr 
times cale as heavie r objects sink inw ards (IChaname & Gouldl 
120021) . Recently, iLiu et al.l (120121) examined the expected 
timing accuracy of pulsars in the GC accounting for radiome- 
ter noise, pulse phase jitter, and the interstellar scintillation 
of the ISM. They found that the 1 hr timing accuracy of SKA 
is expected to be between 10-100 /is for regular pulsars. Our 
results indicate that the necessary accuracy to detect timing 
variations associated to individual 10M Q BHs within 1 mpc 
requires much higher timing accuracy, which might be pro- 
hibitively difficult even with MSPs with a factor of 100-1000 
better timing accuracy. However, the net variations caused 
by a population of these objects is detectable between 2-5 



mpc at these accuracy levels. Remarkably, a 10-100 ^s tim- 
ing accuracy is sufficient to individually resolve or rule out 
the existence of 10 3 M Q IMBHs within 5 mpc from SgrA*. 

As the GW spectrum is rich in strong spiky features at high 
frequencies, it may be possible to isolate GW induced tim- 
ing residuals from other systematic effects. One such effect 
is if the pulsar itself is a part of a binary system. Fortunately, 
however, binaries with orbital periods of years, matching the 
GW foreground of the GC, are very soft, and are not long 
lived near the GC, they are easily disrupted by three-body en- 
counters. Indeed, due to the high velocity dispersion in the 
GC, a ~ 200km/s at D ~ pc, stellar mass binaries are soft 
for orbital frequencies / < (27r)- 1 cr 3 /(GM Q ) - 9,600nHz, 
or orbital period f > 1.2 day. The evaporation timescale 
on which a series of more distant encounters gradually in- 
creases the binary separation, is t e ~ 0.06(j/(G/9alnA) ~ 
2/3 

40/ g Myr, where A is the Coulomb logarithm, and the ion- 
ization timescale to disrupt the binary by a close three body 

encounter is f, ~ 0.04g/(G pa) ~ 70/ g 2/3 Myr (see § 7.5.7. in 
iBinnev & Trema ine 2008). Here we have expressed the bi- 
nary semimajor axis a with the orbital frequency, which is of 
order 0.1 < f% < 10 to match the GW signal. Thus, MSPs 
which form in short period binaries, like in typical LMXBs, 
and become wide, may be expected to typically become sin- 
gle in the GC. However, a few pulsar-BH binaries may form 
through three b ody exchange interactions and m ay be longer 
lived in the GC dFaucher-Giguere & Loebll201 lb . Ultimately, 
multiple pulsars would be necessary to rule out systematic ef- 
fects. 

GW observations with pulsar timing in the GC could be 
combined with other observable channels to map the GC. An 
IMBH orbiting around the SMBH would be observable with 
future millimeter VLBI imaging with EHT (Bro derick et al] 
201 li). Conversely, if those observations reveal an IMBH, the 
inferred orbital parameters could help in identifying the GW 
counterpart with pulsar timing. 
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APPENDIX 
ECCENTRIC SOURCES 

Here we present the mathematical derivation of the GW foreground of eccentric sources in the GC. We start by reviewing the 
eccentric waveforms and GW spectra, then calculate the GW background of an eccentric population of periodic sources and burst 
sources, respectively, and finally discuss simple estimates of the signal-to-noise ratio of timing residuals for individual eccentric 
sources. 



Waveform 

An individual eccentric source with semimajor axis a and eccentricity <?, generates a GW strain composed of discrete upper 
harmonics with frequency /„ = nf 0I \,, 

oo 

h(a,e,t) = ^/i„(a,e;/„)e 27r ' 7 "' (Al) 
«=i 

where 

2 , 

h n (a,e;f,) = -\J g(n,e)h (a), (A2) 
n 

where ho(a) = y / 32/5M.TO, t /(Z)a) is the GW strain amplitude for circular orbits, and 



g(n,e) = — 



J n - 2 - 2eJ n -i + ^J„ + 2eJ n+i -J n+ tj + ( 1 - e 2 )(J„- 2 - 2J„ + J„ +2 ) 2 + 



(A3) 



Here J, = Jj(x) is the ith Bessel function evaluated at x = ne dPeters & Mathews! 19631) . and we have RMS averaged over the binary 
inclination. The Fourier transform of h(a,e,t) measured for some time T is then 

ka,e,f) = Y, h n{a,e-J n )Tw{f,f n ,T), where w(f,f n ,T) = ^ffi^F^ ■ < A4 > 

^ 27r(/-/„)r 

Here w(f,f„,T) is the Fourier transform of a window function of width T and a unit integral. It approaches unity at / = /„ and 
cuts off for \ f-f n \ > 1 /2T, The spectral width of each GW harmonic is A/ ~ 1 /T. 

For circular sources, g{n, 0) = 1 for n = 2 and otherwise. For eccentric sources, the dominant frequency harmonic is f v = ttp/ OT b, 
the inverse pericenter passage timescale, 



rip(e) = ceil 



1.15 (1+e)1/2 



(l + e) 1/2 

so that / p (fl,e)~ ^— ^/ or b(a), (A5) 



(l-«)3/2 

where ceil(x) is the nearest integer larger than x dO'Learv et alj |2009). The emitted GW spectrum is broadband with a maximum 
near / p , where 90% and 99% of the power is between 0.2 / p < / < 3/ p and between 0.1 / p < / < 5/ p , respectively. The g(n,e) 
harmonic weights in Eq. ( IA3b have a maximum near n p (e). 

The definition of the inclination-averaged GW strain amplitude (lA2b can be made more lucid by recalling the definition of the 
GW flux, S = h(t) 2 /16n, and verify if the total power output is consistent with Eq. (16) of iPeters & Mathewsl (11963b . Indeed, 



P = 4nD 2 S = \d 2 1? = ^D^f^hl = A^fl^e) = —^F^ (A6) 
where we have used Kepler's law ui 2 lh = 47r 2 / ( 2 rb = M, / a 3 , the definition of Iiq, 

F(e) = J2g(n,e)= ™ J/2 , and F x (e) = 1 + (73/24) e 2 +(37/39) e 4 . (A7) 

GW background of periodic sources 

Let us estimate the net contribution of many sources to the GW background if observed for time T. For a source with semimajor 
axis a and eccentricity e, the counts in each frequency bin / are given by 

hUa,e;f) = P(a,e;f)f Af = f)/T^«) ^ g (n,e) x { J ^Lis^ ~ ^ "* * > ' ^ 
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where f„(a) = nf mb (a) = n(27rM.) 1 (a/M,) 3 / 2 is the GW frequency of the « th upper harmonic for a fixed semimajor axis. Con- 
versely, the range of semimajor axis, a„±0.5Aa n , for which the « th harmonic contributes to the frequency bin between /± 0.5 Af, 
is 

2 /3 

a„=M.(2irM.f/nT 2 / 3 = fl 2 . (A9) 

where Aa„ = \da„/df\ Af = |a„/(/T). 

Now let us assume a phase space distribution, in which the number of objects in the neighborhood of (a,e) is d 2 N = 
(d 2 N / dade)dade . The number of sources that contribute to the n harmonic is AN,, = f de(d 2 N / da„de) Aa„. The GW sig- 
nal of all sources is then 

(A10) 

a=a„ 

For periodic sources, « max is set by the condition that objects are observed for at least one orbit, 1 /f OI b(a n ) = n/f< 1 /T, implying 
that n max = fT, We shall consider the contribution of burst sources, on larger radius orbits observed for only a fraction of the 
orbit, separately below. Rearranging and using Eq. ( |A8l l, (d 2 N / dade) = A-Ka 2 n+{a)(p{e), and da/df = \aj f, gives 

^LerC/) = yZ~r a n n*(«n) h o( a n) / de (p(e) ^g(n,e) . (All) 
tT 3 Jo n 2 

Let us assume that oc a ~ a . T hen using /zo(fl) oc a -1 and Eq. ( lA9l l. we find that the RHS is proportional to a)~ a . Now let us 
express a„ with «2 using Eq. ( IA91 . 

^c,per(/) = yfl2«*(fl2)/lo(fl2)^per(/) (A12) 

where 

/T /7\ (2+Q)2/3 

*per(/) = J Q deYl [-) vWgfae) ■ (A13) 

Thus K pel (f) = at / < 1 /T, and increases mono tonic ally, and asymptotes a constant for n max = /T — > oo. One can show^that 
this constant is insensitive to the highest eccentricity sources if a > 1/2, which is expected to be sa tisfied . For a = 2 and a 
thermal distribution of eccentricities (pie) = 2e, we get K pa - ~ 0.8 for fT 3> 1. This together with Eq. ( IA12b shows that the net 
GW spectrum of a continuous population of eccentric sources is very similar to that of circular sources. However, the signal is 
much different for individually resolvable sources, as they are comprised of many upper harmonics 



',perv 



(/)= / deJ2 



d l N 
dade 



da 



Af)hlAa,e;f) 



GW background of burst sources 

GW bursts are generated by objects which make only one close approach near the SgrA* during the observation. These sources 
are on eccentric orbits with orbital time f~^ h exceeding T, but for which the pericenter timescale is less than T and the orbital 
phase is such that pericenter passage occurs within the observation. The later condition means that only a l/(J p T) fraction of 
all such sources will contribute in the observation time. For a fixed measurement frequency, /, therefore / 01 b = f/n < l/T, 
frip/n > 1 /T. Therefore fT <n< n p fT and the fraction among these sources that contribute is n/(n p fT). 

Repeating the derivation for the net GW burst background over time T, Eq. (lA8J4A10l i. we get 



where 



hlbif) = -^aln*(a 2 )h 2 {a 2 )K b (f) (A14) 



rl "pf T n / 2 \(2+<*)2/3 

K b (f)= deV— - - ip(e)g(n,e). (A15) 

Je mm ^ T n pf T \ n > 



Note that resolving the burst source also requires that the sampling frequency / max = 2/ At and timescale between observations 
Af to satisfy f p < / max and At/2 < f~ l . This implies that > jn p fAt has to hold. This requirement, however, is typically 
already satisfied for most sources if Af > 1 wee k, the number of sources on so eccentric orbits with such small r p , is close to zero. 
Comparing K pel {f) and K b {f), Eqs. ( IA1 3b and JA15t . shows that the net GW signal of periodic sources exceeds the contribution 
of GW burst sources. 



5 To see this, consider an approximate sharply peaked signal around / p , 
for which g(n, e) ~ S n<n ( e )F(e), where 8jj = 1 if i = j and otherwise. 
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Signal to noise ratio 

Here we present simple estimates on the scaling of the GW signal and timing SNR for eccentric sources using elementary 
functions. This is useful not only because the exact signal presented above is algebraically complicated, but also because the total 
SNR for individual sources is not represented well by h c (f) in a single frequency bin, but it is sensitive to the coherent sum over 
many orbital harmonics. 

Let us define an effective GW strain amplitude using that the GW power at peric enter passage, P p , is the average power times 
the fraction of time the source spends near pericenter passage, i.e. n p , given by Eq. (IA5t . Thus, 

P p = jDW p h 2 p = n p P= ^^(l + e^F^e). (A16) 

p 

This amounts to the root-s um-s quare of individual frequency harmonics in the spectrum for single sources, without averaging 
over the full orbit (c.f. Eq. (IA6I 1). From this we get, that the effective strain at close passage is 

(>yP)'/ 2 [32 M.m* 1/2 

hp = 2— jr = \ — F l (e) (A17) 

Dlo p V 5 Dr p 

where uj p = 2irf p and we have used uj p =M 1//2 (1 + e) l l 2 r^~ . The pulsars are sensitive to the time integral of the strain (see § 14.21 ), 

K,r ^1/2 /32 lJf| / V Ti;2 (l- f )'/ 2 rl/2 , 



p 



n c = ^(f mb T)^ = \~^T^)^Ft%). (A18) 



5 7T Da 1 / 4 (1+e) 1 / 2 1 



The SNR is then S/N = H c / 6 to where 6 to is the timing noise over time / 1 . Note that the RMS of the eccentricity dependent terms 
is 0.4 38 in Eq. ( IA18I ), if the eccentricity is drawn from a thermal distribution tp(e) = 2e between < e < 1 dBinnev & Tremainel 
120081) . The SNR of the timing residual is not very sensitive to the semimajor axis, or the maximum eccentricity in the cluster. 
For fix ed a, eccentric sources contribute to the net timing residuals at a similar level as circular sources. 
Eq. dAl 8b is only applicable for periodic sources, i.e. if f or bT > 1. For a single GW burst source with l// p < T < 1 //orb, 

_ h p _ 1 [32 M l J 2 m, rf 1/2 
oj p 2n V 5 D {l + eyl 1 

provided that the orbital phase coincides with pericenter passage during the observation. This shows that for burst sources with 
fixed pericenter distance, the total timing residual is weakly dependent on eccentricity, for fixed pericenter distance. The number 
of sources is constant as a function of r p for fixed e, but the fraction of sources that are near pericenter passage at any given instant 
decreases quickly for smaller r p proportional to l/n p . Therefore the net contribution of burst sources scales as 

u \mu u 1 tW2 a h P, -ml/2 a 1 M? 4 m* T x l 2 n rp Ar p 1/2 

H c ,b = ANHcAb = HcAb(u p T) 1 n rp Ar p = —{u p Ty'-n rp Ar p =-\ — rn- — ; TTF> F i ^ ( A2 °) 

u)„ 7r V 5 r>rJ (1 + eyi L 



The spectral density within a frequency bin follows from Ar p = (2/3)uj p ^ Aui p implying that % c u cx r p l ^ 4 Ar p oc u^ 2 Auj p . The 
average timing spectral density of burst sources decreases quickly toward higher frequencies. 
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